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A  PROBABILISTIC  THEORY  OF  ANTI-SUBMARINE  WARFARE  MODELS 
DEVELOPED  IN  TERMS  OF  CONGESTION  THEORY 


By 


Brian  W.  Conolly 


>  ABSTRACT 
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This  report,  which  is  methodological,  develops  a  probabilistic 
theory  that  has  direct  application  to  both  antisubmarine 
warfare  and  congestion  models.  The  theory  is  expressed  in 
congestion  terminology  because  of  the  presumed  wider  knowledge 
and  appeal  of  this  field.  This  results  in  a  simplified 
presentation  of  the  general  theory  of  infinite  service  L 

facility  systems  with  specific  application  to  M./Y/®  and 
X/M/«,some  of  which  have  already  been  studied  by  Takacs  and 
Khintchine.  A  new  result  is  given  for  the  output  of  the 
latter  process.  The  analogy  between  certain  infinite 
service  facility  systems  and  a  single-server  system  with 
queue  length  dependent  service  is  exploited  to  provide  results 
for  the  latter  process,  A  further  new  result  for  the  busy 
period  of  such  a  process  is  quoted.  The  antisubmarine 
applications  are  to  the  formally  similar  models  of  the  number 
of  units  present  in  a  geographical  area,  and  to  the  attrition 
of  an  enemy  submarine  force  subjected  to  a  steady  threat  from 
an  antisubmarine  barrier  that  geographical  or  other 
constraints  compel  it  to  transit, 

a\ 

'  \ 
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INTRODUCTION 


This  report  is  concerned  with  the  theory  of  a  probabilistic 
model  that  may  be  used  to  describe  generalized  situations 
that  arise  in  a  military  context,  in  particular  in  anti¬ 
submarine  warfare.  One  situation  is  as  follows:  Units 
arrive  at  the  boundary  of  a  geographical  area  in  which 
they  propose  to  spend  some  time.  They  may  be  ships  about 
to  cross  an  ocean  or  strait,  singly  or  in  convoy;  or  they 
may  be  submarines  going  on  patrol.  In  a  non-military 
situation  they  may  be  cars  entering  a  car  park,  or 
prospective  clients  arriving  at  a  supermarket.  As  long  as 
the  units  are  not  prevented  from  achieving  this  objective 
they  may  be  described  in  the  language  of  Congestion  or 
Queueing  Theory  as  "customers”  entering  a  service  facility 
with  infinite  capacity,  that  is  to  say,  one  that  can  provide 
service  to  any  unit  demanding  it. 

In  these  situations  one  needs  to  make  statements  about  the 
total  number  of  units  in  the  area  at  any  given  time,  and 
about  the  number  of  units  that  leave  the  area  in  a  given 
time  interval.  These  are  respectively  the  content  of  the 
area,  and  the  output  from  the  area.  In  a  defensive  situation 
one  might  be  interested  in  the  number  of  merchant  ship-  in 
an  area  subjected  to  a  submarine,  or  other,  threat  simply 
because  one  wants  to  know  the  scale  of  protection  it  is 
required  to  provide.  From  the  same  point  of  view,  one  is 
likely  to  want  to  know  the  magnitude  of  an  enemy  threat  in 
the  area  in  order  to  be  able  co  assess  the  level  of  forces 
required  to  subdue  it..  The  output  from  an  area  is  equally 
of  interest,  particularly  if  the  "service  time1'  within  the 
area  includes  the  possibility  of  destruction  as  well  as 
safe  t  rans i t . 
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Numerical  statements  about  content  and  output  are  most  easily 
made  in  deterministic  terms.  Thus,  suppose  the  content  of 
an  area  at  time  t  is  £(t)  units,  while  the  rate  of 
arrival  at  the  boundary  is  X  units  per  day  and  the  mean 
time  required  by  each  unit  to  transit  the  area  is  u  1  days. 
Then  the  mean  number  of  units  leaving  the  area  per  day  will 
be  and  we  shall  have 


dt 


x  -  u5  , 


[Eq.  l] 


so  that  if  X  and  U  are  time  independent  and  §(0)  —  % 0 , 


§(t)  =  50  e  ut  +  P(1  -  e  ut),  [Eq.  2] 

where  p  =  x/u*  The  output  rj(T)  from  the  area  in  the 
interval  (0,  T)  is  then 


n(T) 


§(t)  dt 
o 


=  §0(  1  -  e-uT)  +  p  f(T) 


>q-  3] 


where 


f(T)  =  |iT  -  (  1  -  e  uT)  - 


Predictions  made  on  the  basis  of  such  a  theory  can  be  useful 
as  indicators,  but  what  is  realty  needed  in  a  more  operationally 


realistic  treatment  is  a  description  of  §(t)  and  r|(t) 
in  probabilistic  terms.  This  can  be  achieved  by  studying 
the  basic  quantities 


Pn(t)  =  Pr  §(t)  =  n] 


and 


v 

n 


(t) 


Pr 


n(  t) 


n 


Kith  suitable  initial  conditions.  These  are  not,  however,  the 
only  probabilities  of  interest  in  an  operational  study.  The 
duration  of  time,  for  example,  when  ?(t)  >N,  where  N  is  some 
prescribed  number,  is  also  of  interest,  both  in  relation  to 
targets  and  to  the  threat. 

The  same  model  car  be  used  to  describe  the  losses  to  an  enemy 
force  of  submarines  that,  in  order  to  operate,  has  to  cross  an 
antisubmarine  barrier  where  a  constant  kill  rate  per  transit 
can  be  applied.  Let  ?(t)  be  the  number  of  enemy  submarines 
on  operational  strength  at  time  t,  and  \  the  building 
rate.  If  ail  submarines  are  used,  then  the  rate  of  arrival 
at  the  barrier  is  proportional  to  2, ( t )  ,  and  the  number 
destroyed  per  unit  time  can  be  written  u?(t)-  Thus, 
deterministical ly, 


=  \  -  U?  , 


which  Ls  the  same  formal  model  as  that  given  in  Eq.  1  and 
the  same  probabilistic  theory  will  apply. 
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It:  will  be  noticed  that  the  deterministic  differential  equation 
given  above  is  suggestive  of  radioactive  decay..  Familial  i tv 
with  warfare  models  of'  the  Lanchester  type  and  with  epidemic 
and  struggle-f or -survival  models  in  ecology  will  lead  to  the 
recognition  that  the  model  belongs  also  to  this  family,  though 
with  less  complexity  in  that  we  are  not  dealing  directly 
with  the  interaction  of  two  species.  More  generally  the 
model  is  a  particular  kind  of  birth  and  death  process, 

in  the  language  of  Queueing  Theory  the  problem  of  making  a 
probabilistic  description  of  the  content-  c(0  of  a 
geographical  area  is  known  as  Erlang's  Problem  tor  an 
Infinite  Collection,  Erlang  was  sped  f i rally  concerned  with 
the  application  to  telephone  exchanges  and  in  this  context  ft 
is  the  number  of  calls  in  progress  at  time  t ,  Since  Queueing 
Theory  and  its  terminology  are  widely  known  it  seems  appropriate 
to  proceed  with  the  development  of  the  theory'  in  the  queueing 
theoretic  context,  not  forgetting  the  wide  range  of  applications 
in  other  spheres  of*  activity,  In  the*  terminology  of  queueing 
theory,  ’’units"  are  ref  ei  red  to  as  clients"  or  "customers'. 

Thei  r  "service  time"  is  the  time  spent  in  the  service  facility'' 
or  geographic  area.  The  "input"  to  the  "system"  consists  of  a 
description  of  the  intervals  of  time  separating  the  arrivals 
of  successive  cust  omers/'un  i  t  s  at  the  service  facility.  A 
principal  objective  of  tills  report  ;s  a  simplified  expos  i.t  ion 
<rf  the  probabilistic  theory  of  iniini.tr  capacity  service 
systems  that  are  formally  identical  to  the  military  model 
discussed  in  the  previous  paiagtaphs,. 

A  mathematical  dose  1 i pt  ion  of  seiviee  systems  requites  a 
specification  in  probabilistic  trims,  whore  appi opr i a t e ,  of 
the1  input  ,  the  service  t  i  me  ,  in  some  cast's  the  "queueing 
discipline'"  (the  oidei  in  wii  i  <  li  customers  are  served,  lor 
example'),  as  well  as  of  i  he  numhei  of  service'  points.  On  t  he 
basis  of  these  one  makes  prohabi list  ic  statements  about  how 
long  a  t  pica  I  customer  has  to  wait.  I  lie  length  of  the 
uninterrupted  periods  during  wli  i  <ii  t  iie  i  at  i  1  i  t  y  is  kept  busy, 
the'  interval  oi  time  between  successive  depai  lures  t  coin  i  iie 


system,  and  so  on.  To  these  may  be  added  Erlang's  Problem, 
which  is  also  of  particular  interest  to  the  designers  of 
supermarkets  and  doctors'  waiting  rooms,  as  well  as  to 
those  whose  task  it  may  be  to  protect  shipping  against 
enemy  submarines  about  which  the  only  thing  known  is  a 
description  of  the  input  and  cycle  time  in  probabilistic 
terms . 

Erlang's  Problem  for  a  finite  collection  is  rather  intractable, 
particularly  where  so-called  exact  finite  time  results  are 
concerned.  These  are  required  when  one  wants  to  make 
statements  about  §(t),  which  it  will  be  recalled  is  the 
total  number  of  units  present  at  time  t  after  the  system 
was  initiated,  rather  than  about  how  lj(t)  behaves  in  the 
long  run.  It  turns  out,  on  the  other  hand,  perhaps  at  first 
sight  surprisingly,  that  finite  time  results  are  available 
for  certain  systems  that  offer  infinite  service  potential 
and  that  are,  at  the  same  time,  reasonable  models  for 
practical  information  about  the  upper  limit  that  the  system 
can  achieve.  It  is  therefore  of  obvious  practical  value. 

Before  proceeding  further  it  is  advantageous  to  introduce 
a  fundamental  notation  of  Queueing  Theory ?  The  terminology 
X/Y/N  is  used  to  denote  a  service  system  with  N  servers. 

The  units  requiring  service  are  supposed  to  arrive  at  time 
intervals  whose  distribution  is  denoted  by  X.  If  the 
service  facility  is  full,  they  wait.  When  they  reach  a 
service  point  the  time  required  to  complete  service  has  a 
distribution  specified  by  Y.  The  most  commonly  encountered 
type  of  X  and  Y  is  the  negative  exponential  distribution. 
This  is  always  denoted  by  the  letter  M.  Another  type  of 
frequent  occurrence  relates  to  the  deterministic  case, 
and  is  always  denoted  by  D.  Thus  M/D/ 1  denotes  a  single 
server  system  in  which  arrivals  are  separated  by  time  intervals 
having  a  negative  exponential  distribution,  and  each  customer 
receives  a  fixed  service  time.  In  this  report  we  are  dealing 
with  systems  that  may  be  denoted  by  X/Y/»  . 
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Explicit  time-dependent  solutions  of  Erlang's  Problem  for 
M/Y/oe  and  for  a  special  X/Y/cb  in  which  X  is  governed  by 
a  time  dependent  parameter,  were  given  by  Khintchine  [Ref,  2], 
Subsequently  Takacs  [Ref,  3]  gave  a  solution  for  X/M/a# 
where  arrivals  occur  in  a  renewal  process  with  arbitrary 
inter-arrival  distribution.  The  Takacs  solution  is  expressed 
as  a  Laplace  transform.  Subsequently  it  was  discovered 
[Mirasol  (Ref .  4)  ]  that;  the  output  of  the  M/Y/oo  system  is 
Poisson.  This  aroused  considerable  interest  since  an  infinite 
service  system  can  be  used  to  model  delays  in  arrival  at  a 
service  point  [see  Kendall  (Ref,  5),  for  example].  That  is  to 
say,  one  may  imagine  that  when  a  customer  arrives  he  is  put 
through  an  imaginary  infinite  server  system  before  going  to 
the  serving  point,  and  this  has  the  effect  of  delaying  him. 

The  output  of  the  infinite  server  system  then  becomes  the 
input  of  the  real  service  facility,  and  since  queues  with 
Poisson  input  (essentially  of  M/Y  type)  are  well  understood, 
Mirasol' s  discovery  held  out  some  promise  of  throwing  light 
on  the  notoriously  difficult  delay  problem.  More  recently 
Vere  Jones  [Ref,  6]  has  shown  that  the  output  of  an  infinite 
sequence  of  X/Y/oo  systems,  each  feeding  the  other,  is  Poisson. 

With  the  exception  of  a  hint  in  Kendall  [Ref.  5.]}  no  investigator 
seems  to  have  exploited  the  fact  that  under  suitable  renditions 
the  customers  entering  an  infinite  service  facility  may  be 
thought  of  as  flowing  through  it  independently.  This  fact 
can  be  used  to  simplify  considerably  the  development  of 
the  theory.  It  is  the  purpose  of  this  report  to  demonstrate 
how  this  can  be  done,  and  to  make  an  application  to  the 
waiting  time  problem  for  an  interesting,  incompletely 
investigated  queueing  system  in  which  service  is  geared  to 
demand. 


1. 


THE  METHOD 


We  suppose  that  at  time  t=0  the  infinite  capacity  service 
system  opens  its  doors,  permitting  service  to  begin  on  the  n 
customers  who  are  waiting.  First  we  give  formulae  for: 

(a)  The  probability  a^^t)  that  m  of  the 
original  n  customers  have  been  served 
by  time  t.  This  is  clearly  independent 
of  subsequent  arrivals  who,  because  the 
service  facility  is  infinite,  will  flow 
independently  through  the  system. 


(b)  The  probability  density  function  bnm(t) 

that  the  m^*1  customer  of  the  original  n 
has  been  served  in  the  small  time  interval 


(t,  t+dt). 


Let  the  probability  density  function  (p.d.f.)  of  the  service 
time  s  of  a  customer  be  denoted  by  b(t).  We  make  no 
special  demands  on  b(t)  except  that 


P 00 

j  b(t)  dt  =  1 
o 


and  that  its  moments  exist.  We  denote  the  distribution  function 
(d.f.)  by  B(t),  i.e. 


B(t)  =  [*  b(s)  ds. 
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and  the  tail  of  the  distribution  by  Bc(t),  so  that 


Bc(t)  =  1  -  B(t). 


The  m  customers  to  have  been  served  by  time  t  can  be  selected 
in  ways  and  since  the  probability  that  any  of  the 

original  n  has  been  served  by  time  t  is  B(t),  we  have 
immediately 


“  (m)  [™T  [Bc(t> 


n-m 


[Eq„  4] 


"th 

To  deal  with  b  (t)  we  note  that  the  m  customer  who 

nm 

has  been  served  in  (t,  t+dt)  can  be  chosen  in  n  ways. 
The  p.d.f.  of  his  service  time  is  b(t).  In  addition  we 
require  that  m-1  out  of  the  remaining  n-1  are  served 
before  t.  Thus 


b 

nm 


oo 


-  («:!)  [»<*>]”" 


n-m 


b(t) 


[Eq.  5] 


Writing  D 


for 


_d_ 

dt 


it  is  then  clear  that 


b 

nm 


(t) 


(2) 


n-m 


[Eq.  6] 


Now  we  consider  arrivals  that  occur  subsequently  to  t  =  0. 
Again,  because  the  service  has  infinite  capacity,  these  can  be 
dealt  with  independently  of  the  fate  of  the  original  n 
customers  in  the  system.  We  concern  ourselves  with 

(a)  the  probability  r  (t)  that  during  (0,  t) 
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n  arrivals  take  place,  of  which  m  have  been 
served  before  time  t,  (not  necessarily  in  order 
of  arrival) ; 


and 

(b)  the  associated  probability  density  function  qnm(t) 
thatin(0,t)  n  arrivals  take  place,  of  which 
the  has  been  served  in  the  small  interval 

(t,  t  + dt)  (not  necessarily  in  the  order  of 
arrival) . 

We  denote  the  p.d.f.  of  the  intervals  between  arrivals  by 
a(t),  which,  as  for  b(t),  we  suppose  to  be  a  properly 
behaved  function  with  moments.  We  also  write 


A(t) 


and 


a( s)  ds  . 
o 


Ac(t)  =  l-A(t). 


Given  that  t  =  0  is  an  arrival  instant  (one  might  suppose 
that  the  latest  of  the  customers  waiting  at  t  =  0  arrived 
at  that  moment)  we  clearly  have 


r 


T  t 

a(s)  A  (t-s)  B  (t-s)  ds, 
J 

o 


and 


rn(t) 


■>  t 

a(s)  A  (t-s)  B(t-s)  ds. 

0  c 

o 


[Eq.  7] 


[Eq.  8] 


10 


Adding  these  we  obtain,  as  might  be  expected,  the  probability 
of  exactly  one  arrival  in  the  interval  (0,  t) .  The 
explanation  of  Eqs.  7  and  8  is  as  follows.  We  suppose  that 
the  one  and  only  arrival  occurs  in  the  small  time  interval 
(s,  s+ds).  The  probability  of  this  event  is  a(s)  Ac(t-s)  ds. 
In  the  remaining  part  of  the  interval,  t  -  s,  Eq.  7  requires 
that  the  arrival  is  not  served,  [probability  Bc(t-s)],  or 
in  Eq.  8,  is  served,  before  epoch  t  [probability  B(t-s)]. 
Adding  over  all  possible  s  completes  the  explanation. 
Generalizing  to  r  (l  Sm£n-l)  we  have 


r 

nm 


(t) 


f  a(s)  B(t  -  s)  r  ,  .  (t  -  s)  ds  + 

J  ’  n-1,  m-1  v  ' 


a(s) 


B  ( t  -  s) 


n-1,  m 


(t  -  s)  ds. 


[Eq.  9] 


The  explanation  is  again  based  on  a  first  arrival  in  the 
interval  (s,  s  +  ds).  Either  this  arrival  is  served  before 
t  or  he  is  not.  Invoking  the  independence  of  flow  concept 
we  see  in  the  first  case  that  during  (s,  t),  n-1  arrivals 

must  occur  with  m-1  services,  and  in  the  second,  n-1 
arrivals  and  m  services.  For  m  =  0  and  n  we  obtain 
single  integrals: 


rn0(t) 


1  t 

a(s)  Bc(t  -  s)  rn_lj0^t-S^  ds’ 
o 


[Eq.  10J 


r 

nn 


(t) 


t 

a(s)  B(t-s)  rn_x  n_j(t-s)  ds. 
o 


[Eq.  11] 
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With  the  definitions 


r 


00 


(t)  =  Ac(t) 


and 


0 


for  m  >  n,  and  either  m  or  n  negative, the  validity  of  Eq.  9 
can  be  extended  to  all  m  and  n,  thus  including  Eqs.  7 ,  8, 

10  and  11. 

Turning  to  the  corresponding  densities  q^m(t)  we  use  a 
similar  argument.  This  function  has  no  meaning  for  m  =  0. 

For  m = 1  we  obtain: 


q^Ct)  =  a(s)  [b(t-s)  r^j^  Q  (t-s)  +Bc(t-s)qn_1  t  (t-s)]  ds; 
o 

[Eq.  12] 


qnm(t)  = 
nm 


i(s)  Tb(t-s)  r  ,  (t-s)+B(t-s)q  ,  (t  -  s)  + 

L  n-1,  m-1  ’  cv  ’  n- 1 ,  m  ’ 


+B(t-S)qn  lj|n  l(t-s)]  ds; 


[Eq.  13] 


for  2  £  m  £  n-1  ; 


q  ( t) 
Hnnv  ' 


a(s)  [b(t-s)  rn  n  l  (t-s)+B(t-s)qn_1^n_1(t-s) J  ds. 


[Eq.  14] 
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The  explanation  of  Eqs.  12,  13  and  14  is  similar  to  that  for 
Eqs.  9)  10  and  11.  The  first  arrival  occurs  in  (s,  s  +  ds) 
and  we  then  have  to  consider  whether  or  not  he  is  the 
customer  to  be  served  in  (t,  t  +  dt).  In  the  case  m  =  1 
he  either  completes  in  (t,  t  +  dt)  or  after  t|  in  the 
case  m  =  n  he  either  completes  in  (t,  t  +  dt)  or  before  t. 
For  2  smsn-1  he  may  complete  in  (t,  t  +  dt),  before  t 
or  after  t.  There  are  no  other  possibilities. 

Again  a  suitable  definition  permits  the  validity  of  Eq.  13 
to  be  extended  to  all  m,  n.  This  is 


q 


n,m 


(t)  =  0 


for  m  =  0  and  m>n  . 

It  will  be  noted  that  subject  to  suitable  conditions  these 

formulae  for  the  probabilities  r  and  densities  q  are 

1  nm  nm 

valid  for  arbitrary  arrival  and  service  intervals  of 
recurrent  type. 

One  can  now  turn  to  Erlang’s  Problem:  that  of  writing  down 
the  probability 


p(^(t)  =  Pr[?(t)  =  k  j  5(0)  =  N]  , 


where  it  will  be  recalled  that  ?(t)  (tsO)  is  the  non-negative, 
integer-valued  stochastic  process  that  describes  the  number  of 
units  in  the  system  at  time  t  after  the  process  was  initiated, 
and  N  is  the  number  of  units  present  at  t  =  0.  There  are 
two  ways  of  doing  this  and  it  will  always  be  supposed  that 
at  t  = 0  all  services  begin  simultaneously,  and  an  arrival 
can  be  regarded  as  just  having  occurred. 


i 


I 
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First  consider  the  case  OsksN.  It  is  then  possible  to  achieve 

the  condition  |(t)  =k  without  any  arrivals  during  (0,  t)  . 

The  contribution  to  is  ^us  Ac(t)a^  jj  ^(t)*  When 

arrivals  do  occur  it  is  still  possible  that  the  k  remaining 

belong  to  the  original  N.  The  condition  is  that  all  n 

arrivals  (n  a  l)  are  served  before  t.  The  contribution 

to  P(k'(t)  is  thus  “  _k(t)  r  r  (t). 

’  n^l 

Now  consider  that  all  but  one  of  the  arrivals  in  (0,  t)  have 

been  served  before  t.  Then  k  -  1  of  the  original  customers 

must  remain,  and  the  contribution  to  becomes 

a„  i  .  2?  r  „  (t).  Generalizing,  if  s  customers'  out 
N,  N-k+1  .  n,  n-1  v 

’  nil  ’ 

of  the  new  arrivals  remain  at  t,  then  k-s  originals  must 

remain  too.  In  this  case  we  have  a  contribution  of 

a  w  1,1  (t)  2  r  (t)  while  0  £  s  £  k.  This  argument  also 

W  •  rl“KTS  n «  n-s 

9  n=s  9 

depends  on  the  independent  flow  concept.  We  have  finally 


P(S)(t)-Ac(t)  aN  N_k(t)  aN,N-k+s^)n2Jx(1;s)  rn,n-s<t> 


[Eq.  15] 


for  OsksN.  With  the  definition  r^^Ct)  =  Ac(t)  note  that 
the  first  term  in  Eq.  15  can  be  included  in  the  sum.  Thus 
Eq.  15  becomes 


p(k)(t^  “  S  aN,N-k+s(t)  ^  rn,  n-s(t)* 
s=0  ’  nis  ’ 


[Eq.  15a] 


When  kiN+1,  arrivals  must  take  place.  For  simplicity 
consider  k=  N+l.  If  none  of  the  original  N  has  been 
served  before  t  then  exactly  one  subsequent  arrival  must  be 
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present  at  time 

,(t)  E 


NO' 


nss+l 


n, 


t.  The  contribution 

,  ( t)  .  Thus 
n-s-1  '  ’ 


to 


(n) 

N+l 


(t) 


is  then 


s=0  nSs+1 


?  i  (t)  . 

n,  n-s-1  '  ' 


[Eq.  16] 


It  is  not  difficult  to  extend  the  argument  to  obtain  the  general 
formula 


N 

z 

s=0 


E 

n^s+k 


r 

n, 


n-s-k 


(t), 


[Eq.  17] 


for  k  2  1. 

Noting  from  the  accepted  definition  of  binomial  coefficients 
that 


for  m>n  and  negative  n  or  m,  we  then  have  automatically 


a 


0 


for  such  values  of  m  and  n,  and  thus  the  validity  of  Eq.  15a 
can  be  extended  to  all  k(i  0).  This  provides  the  solution  to 
Erlang's  Problem  in  some  generality.  The  argument,  though 
correct,  is  rather  dull  and  the  formulae  on  their  own  are  not 
very  exciting. 
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A  more  aesthetically  attractive  formulation,  which  emphasizes 
the  independence  of  flow  concept,  leads  to  the  following  set 
of  equations  for  p^^(t); 


PN-k(t>=aNk(t> 


N-k 

A  ( t)  +  Z 
C  n=0 


aN,  ktn 


a(s)  p(i}(t-s)  ds, 


[Eq.  18] 


for  0  «k  sN; 


(N) 

pN+k 


(t) 


N  p  t  (  .  \ 

aN,  n(t)  a(s)  pk+„(t's)  ds- 
n=0  »  ^ 


[Eq.  19] 


for  0  <  k  <  »  . 

We  consider  first  Eq.  18.  Here  the  number  of  units  present  at 
time  t  is  at  most  equal  to  N,  the  number  originally  present. 
There  may  have  been  no  arrivals,  in  which  case  provision  has  to 
be  made  for  the  right  number  of  departures  from  among  the 
initial  N.  This  accounts  for  the  first  term  on  the  rifoht 
hand  side.  The  second  term  deals  with  the  case  where 
arrivals  have  taken  place.  We  now  argue  that  because  the 
service  facility  has  infinite  capacity,  the  first  arrival 
creates  a  situation  in  which  a  second  infinite  service  facility 
is  initiated  with  a  single  unit  present.  We  then  have  to 
provide  that  the  sum  of  the  number  remaining  at  time  t  out 
of  the  initial  batch  of  N,  plus  the  number  remaining  in  the 
second  system  initiated  with  a  single  unit  at  time  s,  has 
the  correct  value. 

When  the  number  present  at  time  t  exceeds  the  initial  number  N, 
arrivals  have  to  occur.  The  explanation  of  the  single  term  on 
the  right  hand  side  of  Eq.  19  is  then  identical  with  that  of 
the  second  term  on  the  right  hand  side  of  Eq.  18, 
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A  number  of  steps  can  be  taken  to  verify  the  correctness  of 
these  general  formulations  of  the  solution  of  Erlang;s  Problem 
for  an  infinite  service  facility.  One  can,  for  instance, 


show  that 


ksO 


’  V  1  f 


as  it  should.,  for  all  Nil.  In 


addition  one  can  demonstrate  the  equivalence  of  Eq.  1.5  with 
Eqs.  17,  18  and  19  by  direct  substitution.  It  would  be 
tedious  to  continue  in  full  generality  and  we  shall  therefore 
proceed  in  the  framework  of  particular  systems. 

fi 

It  seems  worth  pointing  out  that  the  theory  permits  one 
to  handle  units  that  require  different  grades  of  service. 

Suppose  for  example  that  the  input  consists  of  a  variety  of 
units, some  having  inter-arrival  interval  density  function  a1(t), 
others  a3(t),  and  so  on,  and  that  service  on  the  al(t)  stream 
has  density  function  bx(t),  bs(t)  on  the  ag(t)  stream,  etc. 
Provided  that  t  =  0  marks  an  arrival  and  beginning  of  service 
epoch  for  each  stream  of  units,  the  different  types  can  be 
allowed  to  flow  independently  through  the  infinite  system. 

The  number  ®(t)  of  units  present  at  time  t  is  then,  with 
obvious  notation,  given  by 


S(t)  =  ?4(t)  +  §a(t)  + 


and  Pr£§(t)  =  nj  is  the  convolution  of  ihe  individual 
probabilities  given  by  the  theory.  Thus  one  can  with  the  same 
theory  describe,  for  example  a  mixture  of  slow  and  fast  ships 
crossing  an  area. 
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A  number  of  steps  can  be  taken  to  verify  the  correctness  of 

these  general  formulations  of  the  solution  of  Erlang's  Problem 

for  an  infinite  service  facility.  One  can,  for  instance, 

show  that  Z!  =  1,  as  it  should,  for  all  Nj  1.  In 

ksO 

addition  one  can  demonstrate  the  equivalence  of  Eq.  15  with 
Eqs.  17,  18  and  19  by  direct  substitution.  It  would  be 
tedious  to  continue  in  full  generality  and  we  shall  therefore 
proceed  in  the  framework  of  particular  systems. 

It  seems  worth  pointing  out  that  the  theory  permits  one 
to  handle  units  that  require  different  grades  of  service. 

Suppose  for  example  that  the  input  consists  of  a  variety  of 
units, some  having  inter-arrival  interval  density  function  a^(t), 
others  a2(t),  and  so  on,  and  that  service  on  the  a.^t)  stream 
has  density  function  b1(t),  bg(t)  on  the  ag(t)  stream,  etc. 
Provided  that  t  =  0  marks  an  arrival  and  beginning  of  service 
epoch  for  each  stream  of  units,  the  different  types  can  be 
allowed  to  flow  independently  through  the  infinite  system. 

The  number  §(t)  of  units  present  at  time  t  is  then,  with 
obvious  notation,  given  by 


5(f)  —  ( t)  +  ?2(t)  +  ...... 

and  Pr£§(t)  =  nj  is  the  convolution  of  the  individual 
probabilities  given  by  the  theory.  Thus  one  can  with  the  same 
theory  describe,  for  example,  a  mixture  of  slow  and  fast  ships 
crossing  an  area. 
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2.  THE  SYSTEM  M/Y/» 

2. 1  Erlang’s  Problem 

This  is  the  system  for  which  the  solution  of  Erlang’s  Problem 
was  first  furnished  by  Khintchine  [Ref.  2],  whose  results  are 
given  here.  The  arrival  pattern  is  governed  by  a  negative 
exponential  distribution  with  parameter  That  is  to  say 


a(t)  =  Xe'U, 

and 


4(t)  =  1  -  e 

A  (t)  =  e~Xt. 
c  ' 


[Eq.  20] 


[Eq.  21] 


Y  denotes  a  general  service  time  distribution  whose  p.d„f.  and 
distribution  function  will  continue  to  be  denoted  by  b(t),  B(t). 
The  simplest  approach  to  the  probabilities  p^^(t)  is  furnished 
in  this  case  by  the  functions  r  (t)  [c.f.  Eqs„  9>  10  and  ll]„ 
Clearly, 


rio(t)  X 


e-Xs  e  B  (t-s)  ds 

c 

o 


=  X 


B  (t-s)  ds 
cv 

o 


=  X  e"Xt  g(t)  , 


[Eq.  22] 
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where 


g(t)  = 


Similarly, 


rXl  =  x 


B  (s)  ds  . 
cv  ’ 


it 


;~Xs  B(t-s)  ds 


X  e‘Xt  f(t)  , 


where 


pt 

f(t)  =  j  B(s)  ds. 


It  will  be  noticed  that  f(t)  +  g(t)  =  t. 


Then  it  is  easy  to  show  inductively  that 


[Eq.  23] 


[Eq.  24] 


[Eq.  25] 


[Eq.  25a] 


r  (t) 
nm 


xn  e-Xt  fm(t)  gn~m(t) 
mi  (n-m) i 


[Eq.  26] 


and  that 


q  (t)  =  £-  D [fm( t )  1 

nm  mi  (n-m) J  L  J 


where  D  =  -j-£  .  Then  from  Eqs.  15  and  17  we  obtain 


[Eq.  27] 


P(k)(t)  =e_Xg(t)  E 

m-max( 0, k-N) 


(„  ^  )BN-k+m(tXBk-m 

vN-K+m'  v  •  c 


[xg(t)] 


m 


mi 


[Eq.  28] 
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for  all  k  a  0,  a  weighted  combination  of  compound  Poisson  terms 
with  time-dependent  parameter  \g(t) .  If  we  put  N  =  1  we  have, 
omitting  arguments, 


[Eq.  29a] 

[Eq.  29b] 


[Eq.  29c] 


[Eq.  29d] 
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Substituting  Eqs.  29a  in  the  integral  in  Eq.  30a  gives 


f*  e\s  -  \g(s)  B(s)  ds 

o 


f*  e*f<s>  de(s) 

J  XL. 


whence  it  is  clear  that  Eq.  30a  is  satisfied  by  Eq.  29a. 
Similarly,  in  the  general  case.  The  coefficient  of  XB  e 
Eq.  30c  is 


Af  (s) 


B„(  s) 


[lg(s)]k  [xg(s)]' 


k+1 


c'  '  kj  +  ( k  +  1 )  [ 


ds 


-  f 

[xg(s) lk+1 

i  r  d 

Xf(s)  L  ;  Jt 

X  J 

0 

e  <k  +  1  < 

_l 

X 


u 


,Xf(s) 


[Xg(s)]k+1 


=  i  PXf(t) 


( k+1 ) J 

[xg(t)] 

( k+ 1 ) J 


k+1 


Thus  the  right  hand  side  of  Eq.  30c  becomes 


-xg  +  R  ft-xg 


B  e  +  B  e 

C  K  • 


.  .  .  .  (i) 

whxcn  xs 
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2.2  The  Output 


Description  of  the  output  process  from  an  infinite  service 
facility  can  be  quite  general,  as  we  shall  now  briefly 
indicate.  Let  n(t)  be  the  number  output  from  an 
infinite  service  facility  in  the  interval  (0,  t)  and  let 


v(^(t)  =  Pr  [n(t)  =  k  |  ?(0)  =  n]  , 


[Eq.  31] 


on  the  assumption  that  t  =  0  corresponds  both  to  the  beginning 
of  an  arrival  interval  and  to  the  beginning  of  service  of  all 
units.  Then  for  k  =  0  we  have  either 


.(») 


(t)  = 


Ac(t) 


an0(t)+an0(t> 


a(s) 


.(1) 


(t-s)  ds 


[Eq.  32] 


or 


(t)  =  Ac(t)  anQ(t) 


+  an0<t> 


mal 


rm0(t)- 


[Eq.  33] 


For  the  M/y/»  system  with  a(t)  =  \e  Eq.  33  gives 

v‘;)(t)  =  B”(t)  e"Af(t)  . 


[Eq.  34] 


Then 


a(t-s)  v^(s)  ds 


=  \e_Xt  eXs  Bc(s)  e"Xf(s)  ds 


=  \e 


-Xt 


eXg(s)  dg(s) 


=  e"Xt(eXg(t)  -  1)  , 


and  it  is  clear  that  Eq.  32  is  satisfied.  The  argument  leading 
to  Eqs.  32  and  33  is  that  none  of  the  original  n  must  have 
been  served  in  (0,  t),  nor  must  any  of  the  subsequent  arrivals. 
In  the  case  of  formulation  by  Eq.  32  it  is  argued  that  the 
first  arrival  initiates  the  activity  of  a  second  infinite 
service  facility.  The  reader  will  readily  verify  that  a  general 
formulation  on  the  lines  of  Eq.  32  gives: 


,(n) 


(t)  =  Ac(t)  ank(t)  +  an.(t)  ^ 


i(t-s)  v£*|  (s)  ds, 


[Eq.  35a] 


for  ksn;  and 

v(k}(t)  =  T  ani(t)  J  a(t-s)  v£X?  (s)  ds 

i  0  A 


[Eq.  35b] 


for  k  >  n. 

In  the  spirit  of  Eq.  33  we  have 


,(n) 


(t)  =  Ac(t) 


nk 


k 

(t)  +  £ 

i=0 


a  •  ( t) 
m'  ' 


m2max(  1  ,k-i) 


r  1  -(e), 

m,k-iv  '  ’ 


[Eq.  36a] 
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for  k  £  n;  and 


T 


i=0 


a  . 
ni 


(t) 


£ 

mak-i 


r 

m} 


k-i 


[Eq.  36b] 


for  k  >  n. 

For  the  system  M/Y/«  these  reduce  to 


min(k, n) 

£ 

i=0 


a  .(t) 
ni 


‘xfCt)]^1 

(k-i)  J 


[Eq»  37] 


for  all  k  s  0.  Thus,  as  for  the  state  probabilities,  the  number 
output  from  the  system  in  (0,  t)  has  for  distribution  a  weighted 
sum  of  compound  Poisson  terms,  but  with  time  dependent  parameter 
\f(t).  This  is  essentially  the  finding  of  Mirasol  [Ref.  4]* 


I 


3. 


THE  SYSTEM  X/M/«o 


k’ 


f- 


3 • 1  Erlang's  Problem 

Although  the  integral  equation  formulation  for  the  output  and 
state  process  seems  more  elegant  it  is  not  clear  that  it  is  as 
suggestive  of  the  solutions  as  the  alternative  expressed  in  terms 
of  the  functions  r  (t) ,  at  least  so  far  as  the  system  M/Y/co 
is  concerned.  The  reverse  is  the  case  for  the  system  X/M/co 
and,  by  contrast  with  M/Y/»,  it  seems  expedient  to  resort  to 
the  Laplace  transformation.  The  density  function  a(t)  is  now 
general,  whereas  b(t)  has  the  specific  form  ye  . 

The  solution  of  Erlang's  Problem  has  been  published,  with  misprints, 
by  Takacs  [Ref.  3]*  It  requires  much  algebra,  which  it  would 
be  tedious  to  reproduce.  We  content  ourselves  with  a  statement 
of  the  result  for  i£(0)  =  1  together  with  an  indication  of  a  method 
for  its  derivation. 

(  n) 

We  drop  the  superscript  on  p  k^(t),  it  now  being  understood 
that  it  is  unity.  We  also  write 


n  oo 


rrk(z)  - 


e"zt  Pk(t)  dt 


for  the  Laplace  transform  of  Pk(t),  and  further  introduce  the  notation 


(^(z)  =  e-(z  +  kuH  a(t)  dt>  ) 


6k(  z)  =  1  "  <*!<(  z)  > 


fk(z)  = 


"  e-(z  +  k^)t  A  (t)  dt; 


[Eq.  38] 
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and  henceforward  drop  explicit  reference  to  the  argument  z. 
Then,  for  kiO, 


"k=  * 

k  i=0 


/k+i\  ,  xi  eia8  *  * ‘“k-i+i  fk+i 

\  k  j  *  *  1  ’ 


Mi  •  •  •  ®k+i 


[Eq.  39] 


in  which  when  k  = 0,  the  numerator  is  fQ  for  i  =  0  and  fj 
for  i  =  1 .  The  generating  function 


n(x)  =  z  tt.  xk 
kso 


can  be  written 


n(x) 


r  (  ,i  CLia«---»k-i+i  fk+± 

izO  Ml-“8k+i 


[Eq.  40] 


Putting  x  =  1  gives 


n(i) 


l 

z  ’ 


which  implies  that  £  Pi,(t)  =  1,  as  it  should. 

k^O 

The  proof  of  Eq.  39  is  carried  out  most  expeditiously  by 
introducing  a  shift  operator  E  that  is  such  that  it  increases 
the  argument  z  of  the  various  Laplace  transforms  by  an 
amount  Thus 

Ok  =  Ek  a0  ,  [Eq.  41] 
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for  example.  Later  it  will  be  convenient  also  to  use  a  backward 
difference  operator  y  which  is  defined  by 


V  =  1  -  E  .  [Eq.  42] 

We  notice  also  that  the  interpretation  of  expressions  like 
(  1  -  tt0  E)"1  is 

i  +  a0  e  +  a0  E(a0  e)  +  a0  E(a0  e)  (a0  e)  + . 

=i  +  a0  e  +  a0  ax  e2  +  a0  e3+  .....  [Eq.  43] 

Finally  we  observe  the  easily  verified  identity 


[Eq.  44] 


Care  has  to  be  ekercised  in  maintaining  the  proper  order  of  the 
operators,  proceeding  systematically  from  right  to  left. 

Now  from  Eq.  18  with  N  =  1  and  k  =  1  we  have 


,  ,  n  t 

p0(t)  =  (l-e-u  )  A  (t)  + ( l-e"wt)  a(s)  p0(t-s)  ds 

c  (J 


which  leads  upon  Laplace  transformation  to 


TT0 


v  |_f0  +  a0 


[Eq.  45] 
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We  write 


Fk 

G0  =  fo  +  Fo  • 


[Eq.  46] 


Then  from  Eq.  45 


G0  f0  +  ao  ^  G0 

(1  -a0  *  )  G0  -  f0  . 


Replacing  V  by  1  -  E  gives 


(«c  +  a0E)Ge  =  f0. 


or 


> 


i  .e. 


G0 


(‘ 


[Eq.  47] 


Next  we  turn  to  the  formula  for  pj(t).  This  is 


pt(t)  =  e_utAc(t)  +  ( l-e“yt)  J  a(t-s)p4(s)+e"ut  J  a(t-s)p0(s)  ds 
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which  gives 


TTX  =  EG0+  VFj 
or 

F:  =  a0  E  G0  +  «Q  V  Fj  t  [Eq.  48] 


and  hence 


F 


1 


[Eq.  49] 


From  Eq.  47 


]• 

after  using  an  identity  of  the  form  of  Eq.  44*  Thus  symbolically 


[Eq.  50] 


For  k  2  2, 


Pk(t)  -<rut 


■>  t  , 

a(t-s)  pk  j (s)ds  + ( 1-e-^  ) 

0 

[Eq.  51] 


1 1 

a(t-s)  Pk(s)  ds, 
0 
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which  leads  to 


30 


which  confirms  the  special  form  of  Eq.  39  for  k  =  0.  This 
result  can  be  cross-checked  for  the  system  M/M/«b  by  putting 
b(t)  =  jje  ^  in  Eq.  29a,  taking  the  Laplace  transform, and 
then  reproducing  the  above  result  for  the  particular  case 
a( t)  =  Xe“Xt. 


3 • 2  The  Output 

The  output  from  X/M/o#  can  be  dealt  with  in  similar  fashion. 

Let  n(t)  be  the  number  of  units  processed  by  the  system  in  (0,t) 
and  recall  the  definition  Eq.  31  with  n=l,  that 


vk(t)  =  Pr  [rj( t )  =  k  |  ?(0)  =  l]  . 


We  use  $^(z)  for  Ehe  Laplace  transform  of  v^Ct)  and 
subsequently  omit  the  argument.  In  addition  we  write 


Hk  ~  ao 


[Eq.  54] 


We  shall  show  that 


H0  =  (l-ft0  E)-1 


[Eq.  55] 


and  that  for  k  2  1 


[Eq.  56] 


This  can  be  more  explicitly,  but  heavily,  expressed  in  various 
ways.  It  seems  best  to  leave  the  formula  in  symbolic  form,  though 
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we  shall  apply  the  check  of  showing  that  it  leads  to 


£  vk (t)  =  i. 

k*0  K 


From  Eq.  3  5a  we  have,  when  n=l  and  k  =  0, 

vQ(t)  =  e-ut  [Ac( t)  +  a(t-s)  v0(s)  dsl 

o 

which  gives 


*0  =  E(fo  +  «o  *o)>  or 
H„  -  a0  E(f„  +  H0). 


Hence 


(1  -  a0  E)  H0  =  a0  E  f0  =  a0  fj 


This  is  Eq.  55»  For  k=l  we  have 


( t)  =  ( 1-e  ^)A  ( t)+(  1-e  ljt)  a(t-s) v0(s)ds+e“ut  a(t-s)vx( 


or 


H!  =  »0  V  [fO  +  HoJ  +  E(K1>  * 


1  + 


1  -  Ot0E 


a0E 


= 


0  i  -  a „E  '  f  o 
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)  ds 


Hence 


H1  1  -  a0E  ’  a°  7  [l  -  ft0E  f o]  > 


which  is  Eq.  58  with  k  =  l.  In  general 


vk(t)  =  (  a(t-s)  vk  l(s)ds+e_|it  a(t-s)vk(s)  ds, 

0  0 

or 


(1  -a0E)  Hk  =  a0  V  Hk-1 


which  plainly  leads  to  Eq.  56. 

As  a  check  we  have: 


\ 


Hence 


2 

ksO 


O-Q 


since  fQ  =  60/z. 

Recalling  Eq.  54  this  means  that 


2 

k^O 


vk(t) 


J 


as  it  should. 


The  author  believes  that  the  results  in  this  chapter  have  not 
been  published  before. 


4. 


THE  THEORY  OF  A  CONSTANT  ATTRITION  ANTISUBMARINE  BARRIER 


Description  in  Terms  of  a  Single  Server  Queueing  System 
in  which  the  Rate  of  Service  depends  on  the  Number  in  the 
Queue . 

4 . 1  Introductory 

An  antisubmarine  model  to  which  the  theory  given  here  is 
applicable  concerns  the  situation  in  which  all  enemy  submarines 
proceeding  from  base  to  operational  area  are  constrained  by 
geographical  or  other  considerations  to  transit  a  region  in 
which  an  antisubmarine  barrier  is  placed.  It  is  supposed 
that  the  barrier  can  exercise  constant  destructive  power  per 
transit.  This  means  that  if  the  enemy  force  contains  n 
submarines  at  time  t,  the  probability  of  a  submarine  loss 
in  the  small  time  interval  (t,  t+dt)  is  jjndt.  This  factor 
coupled  with  the  constant  probability  differential  \  of 
new  construction  accounts  for  the  fluctuations  in  n.  The 
same  model  can  be  used  to  describe  air  defence  situations 
(raiders  attacked  by  missile  defence). 

A  novel  application  of  this  model  in  the  non-military  field 
of  congestion  theory  will  be  defined  in  the  next  section  and 
the  analysis  will  be  prosecuted  in  terms  of  this  model.  The 
terminology  is  not,  in  fact,  of  great  importance,  for  it  is  the 
formal  similarity  of  the  general  model  to  that:  of  an  infinite 
service  system  that  permits  one  to  borrow  from  that  theory, 
thereby  achieving  a  simplicity  in  development  that  would 
be  difficult  to  obtain  by  direct  methods.  In  Sect.  4*3  we  will 
return  to  a  discussion  of  the  antisubmarine  model. 

Research  in  the  theory  of  congestion  is  at  present  embodied  in 
several  books  and  perhaps  about  1500  research  papers,  the 
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majority  of  which  have  been  published  since  1950.  With  relatively 
few  exceptions  these  studies  have  been  made  of  systems  in  which 
the  service  and  arrival  processes  are  independent.  Such  systems 
suffer  from  the  operational  defect  that  when  the  traffic 
intensity,  measured  by  the  ratio  of  mean  arrival  rate  to  mean 
service  rate,  approaches  unity  the  queue  tends  to  elongate 
and  waiting  times  of  the  clients  become  excessive.  On  the 
other  hand,  at  low  traffic  intensities  the  service  facility 
tends  to  suffer  long  periods  of  idleness.  The  obvious  remedy 
for  this  unsatisfactory  state  of  affairs,  which  seems  to  be 
implemented  in  practice,  is  to  arrange,  as  far  as  possible, 
for  service  to  be  matched  to  demand.  (An  alternative  might  be 
to  match  demand  to  service.)  Yet,  curiously  enough,  an  almost 
negligible  amount  of  effort  has  been  expended  on  assessing  the 
effect  of  this  by  analyzing  appropriate  models.  This  is  all 
the  more  curious  at  the  present  time,  in  that  the  effective 
operation  of  computer  time-sharing  systems  may  be  said  to 
depend  on  the  implementation  of  systems  of  this  nature. 

4 • 2  Fundamental  Model 

A  series  of  papers  by  the  author  and  Hadidi  [Refs.  7  to  10] 
have  expounded  the  theoretical  and  numerical  analysis  of  a 
system  in  which  the  allocation  of  service  time  to  a  giv^n 
customer  is  determined  by  the  interval  of  time  between  the 
epochs  of  arrival  of  that  customer  and  his  predecessor; 
substantial  practical  system  improvements  are  demonstrated. 

The  same  authors  have  in  addition  conducted  research  on 
another  system  in  which  the  instantaneous  rate  of  service 
depends  on  the  number  in  the  queue  [Refs.  11  and  12 j. 

As  we  shall  now  explain,  the  analysis  of  this  latter  model  is 
facilitated  by  analogy  with  an  infinite  capacity  service  facility. 

The  system  under  consideration  can  be  explained  as  follows. 

A  simple  M/M/l  system  has  a  single  server  and  is  characterized 
by  independent  negative  exponential  distributions  of  arrival  and 
.service  with  mean  rates  \  and  y  respectively.  We  suppose 


3  6 


i 

I 


f 


& 

$ 


now  that  the  service  parameter  is  on  when  there  are  n 
in  the  queue  and  in  particular  that 


an  =  (n+  1)  u  . 


This  model  then  supposes  that  the  length  of  time  required 
to  serve  a  given  customer  depends  probabilistically  on  the 
number  of  arrivals  during  his  service  time.  A  large  number 
waiting  is  characterized  by  rapid  service.  in  what  follows 
we  shall  further  suppose  that  arrivals  form  a  Poisson 
stream  with  parameter  X  and  that  there  is  just  one  server. 

Let  |{t)  be  the  number  present  in  the  system  (i.e.  in  the 
queue  and  in  service)  at  time  t,  and  write 

[Eq.  57] 


pn(t)  =  Pr  ?(t)  =  n  ♦ 


By  a  familiar  argument,  which  links  |(t  +  dt)  with  §(t) 
via  the  events  that  can  happen  in  the  small  time  interval 
(t,  t +  dt) ,  we  find  that  pn(t)  is  described  by  the  set  of 
equations 


\ 


P0(t)  +  X  P0(t)  =  U  Pi(t) 

P*(t)  +  (X  +  u)  PT(t)  =  Xp0(t)  + 


2(ip2  ( t)  j 


\ 


[Eq.  58] 


Pn(t)  +  (X  +  nii)  Pn(E)  =  XPn_1(t)+(n+l)  UPn+1(t)  j 


where  dots  denote  differentiation  with  respect  to  time.  It 
will  be  recognized  that  these  equations  are  identical  with  the 
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set  that  describes  Pn(t)  for  the  infinite  service  system 

M/M/ ob.  Thus  pn(t)  has  the  same  value  for  both  systems. 

There  is  no  question  of  multiplicative  constants, since 

in  both  cases  72  p  (t)  =  1.  The  results  given  in  Sect.  2.1 
n  n  _Qt 

for  Erlang's  problem  with  b(t)  =  |je  M  are  therefore 
applicable . 

We  shall  now  show  how  to  deal  with  the  problem  of  the 
waiting  time  of  a  customer  who  arrives  when  the  system  has 
been  in  operation  long  enough  to  have  "settled  down"  to  a 
state  of  statistical  equilibrium.  This  means  that  the  number 
in  the  system  is  described  by  probabilities  p^,  which  are 
equal  to  lim  p^(t)  and  have  the  value 


e  P 


n 


[Eq.  59] 


If  we  define  waiting  time  to  include  the  customer's  own  service 
time  and  specify  a  first-come,  first-served,  queue  discipline, 
then  it  follows  that  a  customer  who  arrives  to  find  n  already 
in  the  system  has  to  wait  during  the  residual  service  time 
of  the  customer  being  served,  plus  the  service  times  of 
the  n-1  in  front  of  him,  plus  his  own  service  time.  The 
fact  that  the  service  time  distribution  is  composed  of 
negative  exponentially  distributed  components  means  that 
residual  service  time  has  the  same  distribution  as  complete 
service  time.  Thus  the  new  customer  has  to  wait  for  a  time 
that  has  the  distribution  of  the  sum  of  n+i  complete  service 
times . 


This  leads  us  to  enquire  into  the  form  of  the  joint  probability 
and  density  Ynm(E)  that,  starting  with  n,  it  takes  time  t 
to  output  m  successive  customers.  The  waiting  time  of  the 
customer  who  arrives  to  find  n  in  the  system  is  then  governed 
by  Yn+j  n+j(t) ,  and  the  steady-state  density  function  h(w) 
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of  the  waiting  time  w  of  a  new  arrival,  however  many  he 
finds,  is 


n 

h(w)  -  c-“  Z  It  vn+1  n+1<»>  •  [Eq.  60] 

n^O 


We  now  proceed  to  find  Vnm(^)  by  invoking  the  fact  that  the 
Cn  =  (n+l)u  system  is  formally  identical  with  M/M/eo.  In 
this  infinite  system  is  bo  be  interpreted  as 

referring  to  the  time  t  required  to  output  m  customers, 
no  matter  in  which  order,  when  there  are  initially  n  present. 
If  it  is  assumed  that  at  t  =  0,  service  begins  on  all 
n  customers  (the  shop  opens),  the  problem  can  be  solved 
for  the  system  M/Y/»  in  a  fashion  similar  to  that  employed 
in  Sect.  2.2.  This  we  shall  do.  We  recall  that  T)(t) 
is  the  number  output  from  M/Y/«  in  (0,  t)  and  that, 
by  Eq.  31, 


v(£}(t)  =  Pr  [n(t)  =  k  I  1(0)  =  n]  . 


[Eq.  61] 


Let  T  be  the  time  to  the  completion  of  service  of  the  m 
nm 

customer,  given  that  ?(0)  =  n.  Then  evidently, 


Pr  [Tnm  >  t]  =  Pr  [n(t)  <  m  j  ^(O)  -  n 


n;-  1 

-r  v("*(t). 

k=0 


Equation  62  determines 


V  (t).  since 
Tnnr 


Pr  Tt  >t 
L  nm  J 


i » 

4 


ds . 


[Eq.  62] 
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up#**! ■vv*';  ?  w $  &  1  *  '■r-W’^*#-' 


From  Eqs.  37  and  4  we  have,  for  ksn, 


:j)(t)  -  e-lf<*)  Z  (“)  B^t)  B^V)  li— ill - 

K  i=0  Xl/  C  (k-i)l 

-  e-Xf(t>  c^e^^B^t)  +  x  B(t)]")  ,  [Eq.  63] 


where  C 


Thus 


(x) 


means  the  coefficient  of  x  in 


v(kn)<t> 


-Xf(t)  P 


2rri 


;\zf(t)  |^Bc(t)  +  z  B(  t)  ]  ■ 


dz 

k+1  ’ 


where  F  is  a  contour  in  the  z  plane  enclosing  z  =  0.  Then 


Pr 


T  >  t 
nm 


l  =  e~xf(t)  r  x*f(t)  [Bc(t>  +  zB(t)]n  ^ 

2iri  (  1  -  z)  zm 


[Eq.  64] 


provided  that  r  excludes  z=  1.  Hence 


ce  n- 1 

Hc(t)  =  I  h(s)ds  "  e'°  r,  T^nr  Pr[Tnn>t] 

°t  n^l 


j-p-Xf  ( t)  p 

2rri 


exp 


Xzf(t)  +  p[Bc(t)+zB(t) 


Bc(t)+z  B(t) 


dz 


z(  1  -  z) 

[Eq.  65] 
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Equation  65  gives  the  tail  of  the  waiting  time  distribution  for 
the  single  server  no  =  (n+l)y"  queueing  system  when  b(4) 
is  replaced  by  pe-*^.  The  result  is  delivered  much  more 
expeditiously  by  utilizing  the  infinite  service  system 
analogy  than  by  any  direct  approach  that  we  have  seen. 

[See  Hadidi  (Ref.  12)  where  Eq.  67  is  obtained  by  complicated 
analytical  methods.] 

Hc(t)  can  be  expressed  explicitly  in  terms  of  modified  Bessel 
functions  of  the  first  kind  by  expressing  Eq.  65  in  terms  of  the 
well-known  contour  integral  representations  of  the  latter.  Thus 
one  obtains 


1 

3p 

Hc(t)  =e-5  11 

B  (t)  I  (2  v^rD  +  F  (■ 

v  2  (  1  )  - 

i)  Vi'2*) 

nS0 

i- 

M. 

<•;  . 

% 

where 

& 


JE*. 

& 


&* 


?  =  xf(t) 

11  =  pBc(t)  . 


[Eq.  66J 


Similarly,  the  density  function  is  given  by 


h(t)  =e-5'r!  <[xB(t)Bc(t)  +  b(t)(l  +  pB(t))]  I0(2  J$r\  )  + 


>-  [pb(t)Bc(t)(^)2  +  XB2(t)(|)2]  Ix(2  , 


[Eq.  67] 


which  may  be  evaluated  either  directly  from  the  contour  integral 
obtained  from  Eq.  65  by  differentiation,  or  by  differentiating 
Eq.  66  and  utilizing  known  Bessel  function  relations. 
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Equation  65  can  be  checked  by  setting  t  =  0.  Hc(0)  should  have 

the  value  unity.  Now  f(0)=0,  B  ( 0)  =  1 ,  B(0)=0.  Hence 


Hc(0) 


e  P  P  o/z  dz 

—  Jr  e  KH T 


=  -  l.  r 

2iri  J  ,e  5  ’ 

r‘ 

by  an  obvious  transformation. 

This  is  clearly  unity.  Equation  66  may  also  be  checked  direccly 
in  the  sai,e  way. 

An  explicit  treatment  of  some  analytical  and  practical  aspects 
of  the  an=(n+l)y  single-server  system  can  be  found  in  Hadidi 
and  Conolly  [Ref.  ll]. 


4 . 3  ASW  Application 

The  interpretation  of  this  result  will  now  be  made  in  the 
context  of  the  antisubmarine  model  referred  to  in  Sect.  4.1, 
in  which  all  enemy  submarines  proceeding  to  and  from  their 
operational  area  are  obliged  to  transit  an  antisubmarine 
barrier  that  has  constant  destructive  power  per  transit. 

The  commander  who  asks  how  long  it  will  take  to  destroy  as 
many  submarines  as  the  enemy  now  possesses  can  be  answered  by 
using  the  fact  that  the  distribution  of  the  time  is  1  -Hc(t). 
This  is  the  interpretation  of  the  single-server  waiting  time 
result.  Note  that  unless  X  is  reduced  to  zero  the  enemy 
is  still  likely  to  possess  submarines  even  when  a  number 
equal  to  his  present  strength  has  been  destroyed. 
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A  cognate  question  is;  "How  long  before  I  first  reduce  the 
enemy's  present  strength  to  zero?".  This  is  of  interest  to 
both  offence  and  defence.  More  generally  it  can  be  framed 
in  terms  of  a  reduction  from  the  present  level  to  some  other 
prescribed  non-zero  level.  Such  questions  are  answered  in  the 
terminology  of  stochastic  processes  by  consideration  of  the 
distributions  of  "first  passage  times".  In  congestion  theory 
the  interval  of  time  from  the  epoch  at  which  one  client  first 
occupies  the  system  until  for  the  first,  time  thereafter  it 
becomes  empty  is  described  as  a  "busy  period",  and  one  says 
that  in  this  time  the  system  has  made  a  first  passage  from 
unity  to  zero.  Let  k(t)  be  the  density  function  of  a  busy 
period  for  a  system  in  which  arrivals  and  services  have 
negative  exponential  characteristics.  Then  the  time  required 
for  a  first  passage  from  N  to  n  (N>n)  has  density 
function  k^~n^(t),  the  bracket  superscript  denoting 
convolution.  A  congestion  model  with  negative  exponential 
arrivals  corresponds  to  an  antisubmarine  barrier  model  with 
time-independent  building  rate. 


The  writer  has  analyzed  the  problem  of  the  busy  period  for 
the  M/Y/«  congestion  system,  and  by  analogy  this  is  applicable 
with  negative  exponential  service  time  to  the  barrier  model  in 
the  sense  mentioned  above.  The  result  may  be  expressed  as 
follows . 


Let 


b(t)  =  pe“^  , 
n  t 

B(t)  = 


-ut 

b(s)  ds, 


f(t)  = 


n  t 


B(s)  ds, 


Let 


,  /  .  \  xf  ( t )  . 

h  ( x ,  t  )  =  e  v  -1 


43 


<wwrw! 


where  x  is  some  real  number,  and  let  h^n^(x,  t)  be  the  n-fold 
convolution  of  h(x,  t)  with  itself  (with  respect  to  t). 


Moreover,  let 


h^(x,  t)  =  I  xmhnjn('t^. 


man 


Then  if  Dn  s  d 
dt 


n 


n 


k(t)  =  “T~‘  r  {~)n  Dn+2  [h(n+1) 


n^O 


This  result  appears  to  be  new.  Its  proof  is  given  in  Appendix  A. 

In  terms  of  the  model  of  units  in  a  geographical  area  the 
functions  Ynm(^)  give  information  on  the  time  required  to 
output  m  units, given  that  n  are  present  at  a  certain 
moment.  If  the  units  are  merchant  ships  carrying  strategic 
material  such  information  seems  logistically  pertinent. 
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APPENDIX  A 


THE  BUSY  PERIOD  FOR  THE  INFINITE  CAPACITY 
SERVICE  SYSTEM  M/G/« 


In  the  terminology  of  Queuing  Theory  a  busy  period  (B.P.)  begins 
when  the  service  facility  ceases  to  be  idle,  and  it  continues 
until,  for  the  first  time  thereafter,  it  becomes  idle  again. 

An  infinite  capacity  system  is  supposed  to  be  supplied  with  an 
infinite  number  of  servers.  Immediately  before  a  B.P.  they 
are  all  idle.  At  the  beginning  of  the  B.P.  an  arrival  occurs 
and  immediately  begins  to  be  served.  The  B.P.  continues  until 
once  again  all  servers  are  idle.  It  is  hardly  necessary  to 
add  that  since  the  number  of  servers  is  supposed  infinite  not 
all  servers  will  be  continuously  occupied  during  a  B.P. 
as  here  defined.  Let  n(T)  be  the  number  of  customers 
served  during  a  B.P.  that  begins  at  time  t  =  0  and  terminates 
at  t  =  T.  We  are  primarily  concerned  in  this  paper  with 
the  joint  probability 


kn(t)dt  =  Prj^t  <  T  <  t+dt ,  n(T)  =  n 


[Eq.  A  ,  1  ] 


The  M/G/oe  system  is  characterized  by  an  arrival  pattern  that 
forms  a  Poisson  stream  with  time-independent  parameter  X • 
Thus  interarrival  intervals  have  p.d.f.  Xe  ^ »  A  service 
period  initiated  by  any  one  of  the  infinite  number  of  servers 
is  supposed  to  have  a  proper  p.d.f.  b(t)  and  a  distribution 
function  denoted  by  B(t),  such  that  B(ee)  --  1.  The  duration 
of  a  service  is  independent  of  the  arrival  stream. 

We  also  introduce  the  function 


f(t) 


t 

B(s)  ds, 
o 


[Eq.  A. 2] 
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and  a  generating  function  n(x,  t)  defined  by 


h( x,  t) 


exp  [  X  xf  ( t )  j  -  1  , 


[Eq.  A. 3] 


where  x  can  be  taken  to  be  a  real  variable  restricted  as 
necessary  to  ensure  convergence.  The  commonly  accepted 
notation  p( t)”q( t )  will  be  used  for  the  convolution  of 
two  functions,  viz. 


p(  t)";!q(  t) 


p  t 

I  p(s) 


q(t-s)  ds, 


and  p^n^(t) 
with  itself, 


will  denote  the  n-fold  convolution  of 
viz . 


p(t) 


P(n)(t)  = 


p(n-»(t) 


P(f  )  > 


with  the  convention  that.  p^^(t)  =  p(t).  Then  convolution  of 
h(x,  t)  with  itself  and  with  respect  to  t  produces  a  higher 
order  generating  function  h^  (x,  t)  ,  which,  when  expanded 
in  ascending  powers  of  x  has  coefficients  that  will  be 
denoted  by  h  (t).  Thus 

w  run  ' 


h(n)U,  t) 


I  h  (t)  xm, 
nm  '  ’ 

m^n 


[Eq.  A. 4a] 


and  we  observe  then  from  Eq.  A. 3  that 

r  / 

[_Xf(t)j  /'mi  - 


[Eq.  A. 4b] 
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In  the  following  development  it  will  simplify  writing  to  replace 
expressions  like  dn  p(t)/dtn  by  Dnp(t).  In  general  when 
there  is  no  ambiguity  we  shall  omit  arguments. 

The  objective  of  the  paper  is  to  prove  the  result 


kn(t)  “  ^  [D8(hln>  -D3<h2n>  + . +<->n+‘  ^ 

[Eq.  A, 5] 

The  proof  will  proceed  on  the  following  lines.  Suppose  that 
customers  leaving  the  infinite  service  system  are  labelled 
in  the  order  of  their  departure  (l,  2,  3j  . n,  ...).  This 
sequence  is  a  permutation  of  their  order  of  arrival.  We 
suppose  that  at  t =  0  a  B.P.  begins  and  consider  the  event  E^ 
that  the  customer  labelled  n  has  been  served  in  (t,  t + dt) 
leaving  the  system  empty.  There  are  two  mutually  exclusive 
ways  in  which  can  occurs  either  it  is  the  termination 

of  the  busy  period  beginning  at  t  -  0;  or  there  may  intervene 
one  or  more  idle  periods  before  t.  We  evaluate  the 
unconditional  probability  of  occurrence  of  E^,  and  the 
probability  of  E^  given  the  intervention  of  one  or  more 
idle  periods;  kn(t)dt  is  then  the  difference  between  these 
quantities . 

We  first  require  two  lemmata. 

Lemma  1,  It  is  supposed  that  the  system  is  empty  at  t  =  0. 

Let  rn(t)  be  the  probability  that  exactly  n  arrivals  take 
place  during  (0,  t)  and  that  all  of  them,  and  no  more,  have 
been  served  before  t.  Then 


r(t)=  e  ^  h,  , 
n  ’  In 


(_Eq.  A, 6] 


Proof .  Consider  first  r  (t).  Let  the  first  and  only  arrival 
occur  at  epoch  u.  He  has  to  have  been  served  before  the 
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termination  of  the  remainder  of  the  interval,  which  has 
length  t  -  u.  Then 


rt(t) 


X 


rt  e-xu  e~\( t-u)  B(t_u)  du 

o 


=  e  Xt  f(t)  =  e  Xt  hj  ' 


[Eq.  A  „  7] 


Now  suppose  that  Eq„  A. 6  holds  for  all  n,  and  consider 
r  (t).  If  again  the  first  arrival  occurs  at  epoch  u, 
then  during  t-u 

(a)  this  customer  must  have  been  served, 

(b)  independently,  n  further  customers  have 
arrived  and  been  served. 

The  fact  that  we  have  an  infinite  system  enables  us  to  make 
independent  arrangements  for  the  first  arrival  and  to  separate 
his  channel  from  the  system  as  long  as  it  is  occupied.  Thus 
at  epoch  u  the  initial  conditions  of  the  lemma  are  reoroduced 
and  we  have 


r 


n+1 


(t) 


p  t 

X  e  B(t-u)  r  (t-u)  du 

tJ 

0 


1 


n+1 
n  i 


e~^U  B(t-u)  f  n(  t  -  u)  du 

o 
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by  Eqs.  A. 6  and  A. 4b 


fn  D ( f )  du 


=  e-^  h 


l,n+l  * 


This  completes  the  proof  of  Eq.  A. 6. 

Lemma  2.  Let  t  =  0  be  the  beginning  of  a  B.P.  Let  T  be 
the  epoch  at  which  the  n  service  was  completed,  not 
necessarily  in  order  of  arrival  and  not  necessarily  without 
the  intervention  of  one  or  more  idle  periods  in  the  service 
facility,  and  with  an  empty  system  immediately  after  Tn> 

Let 


vn(t)  dt  =  Pr( t  <  Tn  <  t  +  dt)  . 


Then 


v 

n 


(t) 


[Eq.  A. 8] 


(Note  that)  v  (t)  is  the  P-d.f.  of  the  event  E^  referred 
to  earlier. 

Proof .  This  is  again  inductive.  To  deal  with  vt(t)  we 
require  no  arrivals  in  (0,  t)  and  the  completion  of  service 


51 


of  the  original  customer  in  (t,  t+dt).  Thus 


Vj(t)  =  e  Xt  b(t) 

=  e_Xt  D 2 ( f ) 

e-Xt 

=  D2  (hi;i)  .  [Eq.  A. 9] 

Now  consider  v  .  If  it  is  the  original  customer  who  has 
been  served  in  (t,  t+dt)  we  have  a  contribution  of  amount 

b(t)  r  (t)  =  D(B)  e"Xt  h. 
v  ’  n  In 
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£y-  e~xt  D(fn  B) 


e“Xt  D(fn  B) 


This  proves  Eq.  A. 8. 

To  prove  the  main  result  [Eq.  A. 5]  we  again  suppose  the 
conditions  of  Lemma  2  and,  recalling  earlier  statements,  observe 
that 

Prl  t  <  T  <t  +  dt1  =  PrftcT  <t  +  dt,  no  idle  periods"! 

L  J  L  before  t 

Pr|t  <  Tn  <  t  +  dt ,  at  least  one  idle  |  ; 

period  before  t 

[Eq.  A. 11] 

i.e.  that 

vn(t)  =  kn(t)  +  [Eq*  A.  12] 

where  refers  to  the  second  term  on  the  right-hand  side 

of  Eq.  A. 11.  Knowing  v  ,  if  we  can  determine  1^,  then  we 
shall  have  obtained  k  .  Suppose  the  first  idle  period  occurs 
after  the  termination  of  the  m  service  (l  sm5n-l)  at 
epoch  t-u.  Then  (0,  t-u)  is  a  B.P.  with  m  services. 

We  then  have  to  arrange  for  a  further  arrival  at  epoch 
t-u  +  u>( 0  <  <*>  <  u)  ,  which  reproduces  the  initial  conditions. 
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During  u-<U,  n-m  services  are  required,  with  or  without  idle 
periods,  the  last  terminating  at  epoch  t.  Hence 


tn(t> 


X 


n- 1 

T 


m=l 


■>  t 

k  ( t  -  u)du 
mv  ' 

o 


e 


-XU) 


V 

n-m 


(u  -  id) 


dui 


n- 1 

=  T  k  (t)* 
_ .  m' 


9 


using  Eq.  A. 8.  Now  using  Eq„  A.  5  for  1  imin-1,  we  have 

[D‘(hl,l*hln-l)  + 

+  DS(hl,2*hln-2)-D4(h2,2,,hl„-2)  + 

+  D‘(hl,3  *hln-3)  -D‘(h2,3  *hln-3)  +I>5  (h3,3  *hln-3) 

+  . 

+  D’(hln-1  *  hl,l)  -D‘(h2n-1  *  hl,l)  +I>5  (h3„-l  *hl.l)  ‘ 


,+(.)nDn+I(h 


n-1  ,n-l 


hi(i> 


=  [D“(h2n>  -D‘<h3n>+D5<h4„>  -  Dn+1  (»„„) 

[Eq.  A. 13] 


54 


I'RTiiWT'iiiW'ffln'Tii — r--nr~r~rr^ - r-Tr-1 — rmnniii|)inri(j>i«i  nm  l 


This,  with  Eq.  A. 12  establishes  Eq.  A. 5,  it  remaining  only 
to  show  from  first  principles  that 


Jfc 


kj  ( t) 


D2(h!,l)> 


[Eq.  A.  14] 


which  is  obvious  since  D2(h  )  =  \b. 

*  9  * 

The  p.d.f.  k(t)  of  a  busy  period  of  length  t  is  given  by 


k(t)  =  T  k  (t) . 
nsl 


[Eq.  A.  1.5] 


Using  Eq.  A. 5  we  obtain 

k(t)  =  £  (-)n  Dn+2  [h(n+l)(t)l  ,  [Eq.  A.  16] 

X  n:>0  L  J 

where  h^n^(t)  =  h^n^(l,  t)  . 

As  a  check  we  now  propose  to  obtain  Eq.  A. 16  by  alternative 
means.  The  event  F  consisting  of  the  beginning  of  a  B.P. 
is  recurrent  for  M/G/m  and  the  probability  of  its  occurrence 
in  (t,  t  +  d  t)  is  \P0(t)dt,  where  p0(t)  is  the  probability 
that  a  system  intiated  at  t  =  0  with  one  member  will  be 
empty  at  time  t.  The  interval  between  two  successive 
occurrences  has  density  function  \k(t)  #  e  X\  and  so 


Xp0(t)  =  Xk*e'Xt  +  (Xk*e-XV[xPo(t)]  , 


[Eq.  A. 17) 


since  F  may  occur  for  the  first  time  in  (t,  t  + dt ) ,  or 
may  occur  at  a  previous  epoch. 
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t 


If  tt0(z)  and  K(z)  are  the  Laplace  transforms  of  p0(t) 
and  k(t),  then  Eq.  A. 17  is  equivalent  to 


(X+z)tt0(z) 

K(z)  -  ttt^t 


This  gives  formally 


K(z)  =  z  (-X)n  (X+z)  TT0n+1(z). 
n^O 

Now  it  is  well  known  [e,g,  Khintchine  (Ref.  2)]  that 

pe(t)  =  e'Xt+Xf  B  =  D[h(t)]  , 

where  h(t)  =  h(l,  t) . 

If  we  write  H(z)  for  the  Laplace  transform  of  h(t) 
immediately  from  Eq.  A. 20  that 


tt0(z)  =  (X+z)  H(z+x)/x 

and  hence  that  Eq.  A. 19  gives 

K(z)  =  X"1  I  (-)n  (X+z)n+2  Hn+I(z+X). 
naO 


[Eq.  A  o  1 8] 


[Eq.  A. 19] 


[Eq.  A. 20] 


it  follows 


ijwfticf*  ViWWJf*  ISf.f 


This  transforms  back  into 


I 


5 

i 

i 


k(t)  -  e-^4  X-1  r  (..)"  D,1+2  rh<n+1>l 

n^O  L  J 


which  is  Eq.  A. 16. 

Thus  we  have  a  check  on  the  purely  probabilistic  arguments  of 
the  earlier  paragraphs.  Because  of  Eq.  A.18  and  the  known 
properties  of  p0(t),  we  also  observe  from  Eq.  A.18  that 


P® 

j  k(t)dt  =  1. 
o 

The  mean  B.P.  length  is  perhaps  most  expeditiously  obtained 
from  the  relation,  applicable  to  M/G/»,  that 


E(I) 

E(  BP)  +E(I) 


lim  p0(t)  =  e 

t-*ao 


[Eq.  A. 21] 


where  E(  )  means  "expected  value  of",  and  I  denotes  an 
idle  period.  This  gives 


E ( BP)  =  |  (ep  -  1)  , 


[Eq.  A. 22] 


where  p =  \/E  (service  interval). 


i 
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